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Abstract. The observations of a strong and extended positron-electron annihilation line emission in the 
Galactic center (GC) region by the Spectrometer on the International Gamma-Ray Astrophysical Laboratory 
(SPI/INTEGRAL) are challenging to the present models of positron sources in the Galaxy. In this paper, we 
study the possibility of pulsar winds in the GC to produce the 511 keV line. We propose that there may exist 
three possible scenarios of pulsar winds as the positron sources: normal pulsars; rapidly spinning strongly magne- 
\^ , tized neutron stars (magnetars) in gamma-ray burst (GRB) progenitors; a population of millisecond pulsars in the 

Galactic center region. These e pairs could be trapped in the region by the magnetic field in the Galactic center, 
and cool through the synchrotron radiation and Coulomb interactions with the medium, becoming non-relativistic 
particles. The cooling timescales are lower than the diffuse timescale of positrons, so low energy positrons could 

in . 

annihilate directly with electrons into 511 keV photons or form positronium before annihilation. We find that 
normal pulsars cannot be a significant contributor to the positron sources. Although magnetars in the GC could 
be potential sources of positrons, their birth rate and birth locations may impose some problems for this scenario. 
|' We believe that the most likely candidate positron sources in the GC may be a population of millisecond pulsars 

consistent with the present observational constraints. Therefore, the e ± pairs from pulsars winds can contribute 
significantly to the positron sources in the Galactic center region. Furthermore, since the diffusion length of 
positrons is short, we predict that the intensity distribution of the annihilation line should follow the distribution 
of millisecond pulsars, which should correlate to the mass distribution in the GC. 
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1. Introduction (3 - 6) x 10 42 e+s -1 . Recently, the analyses of the ob- 
servational data by SPI/INTEGRAL with deep Galactic 

Since the first detection (Johnson & Haymes 1973) and center region exposure (longer than 4 Ms) show that thc 

the subsequent identification (Leventhal et al. 1978) of the gpatial distribution of 511 keV line appears cen tered on thc 

Galactic 511 keV annihilation line, the origins of the galac- Galactic center (bulge com p nent), with no contribution 

tic positrons have become a lively topic of scientific debate. from & digk component (Teegarden et al. 2005; Knodlseder 

With the launch of INTEGRAL gamma-ray observatory gt al 20Q5 . churazov et aL 20fj5) C hurazov et al. (2005) 's 

m 2002, the SPI telescope on board provides a strong con- analysis suggested that the positron inje ction rate is up to 

stramts on the morphology and intensity of the 511 keV 1Q 43 e+s -i within „ ^ and an intrinsic annih ii at ion line 

line emission from the Galactic center (Knodlseder et al. w j d ^ n ^ 2 37 ± 25 keV 
2003; Jean et al. 2003). The data analyses show the line 

emitting source is diffuse, and thc line flux within 5° of The potential positron sources include neutron stars 

the GC amounts to ~ (9.9 ± 4) x 10~ 4 photon cnr 2 s" 1 or black holes (Lingenfelter & Ramaty 1983), 56 Co 0- 

(Knodlseder et al. 2003), corresponding to a luminosity decays in supernova remnants (Lingenfelter & Ramaty 

of - 10 36 erg s -1 . The high line luminosity requires that 1979; Ellison, Frank & Ramaty 1990), other radioac- 

the positron injection rate in the GC should be around tive nuclei formed by nucleosynthesis in supernova, nova, 

red giants, and Wolf-Rayet stars (Ramaty et al. 1979), 
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(Kozlovsky et al.1987), pair production by gamma-ray 
photons interacting with starlight photons in the in- 
terstellar medium (Mastichiadis, Protheroc & Stephen 
1991), and gamma-ray cascades at the polar caps of pul- 
sars (Harding & Ramaty 1987; Boulares 1989), electron- 
positron pairs produced by the pulsar winds (Chi, Cheng 
& Young 1996), and probably gamma-ray bursts (GRB, 
Lingenfelter & Hueter 1984). However, the recent results 
from SPI/INTEGRAL on the 511 kcV annihilation line 
emission in the Galactic center show a diffuse source with a 
high line luminosity, and present a challenge to the present 
models of the origin of the galactic positrons. Hence, Casse 
et al. (2004) argued that normal supernova explosion in 
the GC cannot contribute to the required positrons, but 
hypernovae (Type Ic supernovae/gamma-ray bursts) in 
the Galactic center may be the possible sources of galactic 
positrons. In addition, annihilations of light dark matter 
particles into pairs (Boehm et al. 2004) have been also 
suggested to be the potential origin of the 511 keV line in 
the GC. 

In this paper, we will study the possible contribution to 
the positrons in the Galactic center from pulsar winds. We 
suggest that there are three possibilities of pulsar winds 
as positron sources in the Galactic center: normal pul- 
sars (e.g. Crab and Vela like pulsars); the rapidly spinning 
strongly magnetized neutron stars which are possible GRB 
progenitors (Usov 1992); a millisecond pulsar population 
(Wang, Jiang & Cheng 2005). In §2, we will first study the 
possibility that electron-positron pairs can be produced by 
the pulsar winds. In §3, we discuss electron-positrons pro- 
duced in three possible scenarios in detail. These pairs 
could be trapped in the Galactic center by the magnetic 
field. In §4, we find that the diffuse timescale of electron- 
positrons in the GC is much longer than the synchrotron 
cooling timescale and energy loss timescale by Coulomb 
interactions in the medium, so that pairs will become 
non-relativistic, and then annihilate into 511 keV emission 
lines. The summary and discussions will be shown in §5. 

2. Electron-positron winds from pulsars 

It is well known that relativistic particles from pulsar 
winds interacting with the interstellar medium form the 
synchrotron wind nebulae (see the reviews by Arons 1998; 
Kaspi, Roberts & Harding 2004). Chi, Cheng & Young 
(1996) have proposed that the bulk of the comic positrons 
could originate from pulsar winds. In this paper, we will 
consider electron-positron pair production which occurs 
in the pulsar outer-magnetospheric region, the so-called 
outer gap (Cheng, Ho & Ruderman 1986). 

Zhang & Cheng (1997) have developed a self-consistent 
mechanism to describe the high energy radiation from 
spin-powered pulsars. In the model, relativistic charged 
particles from a thick outer magnctosphcric accelerator 
(outer gap) radiate through the synchro-curvature radia- 
tion mechanism (Cheng & Zhang 1996) rather than the 
synchrotron and curvature mechanisms in general, pro- 
ducing non-thermal photons from the primary pairs 



along the curved magnetic field lines in the outer gap. 
The criterion for the existence of the outer gap requires 
the fractional size of the outer gap / < 1, which is the 
ratio between the mean vertical separation of the outer 
gap boundaries in the plane of the rotation axis and the 
magnetic axis to the light cylinder radius. / is limited by 
the pair production between the soft thermal X-rays from 
the neutron star surface and the high energy gamma-ray 
photons emitted from the outer gap region, and can be 
approximated as / ~ 5.5P 26 / 21 i? 12 4 ^ 7 (Zhang & Cheng 
1997), where P is the spin period, and -B12 is the surface 
magnetic field in units of 10 12 G. 

The pair production mechanism is a synchrotron pho- 
ton cascade in a strong magnetic field. According to 
Halpern & Ruderman (1993), each primary electron- 
positron from an outer gap will have an energy 



E p = j p m e c 2 = 5.7 x 10 12 P 1/3 eV 



(1) 



before it strikes the neutron star surface, where 7 P is the 
Lorentz factor of the primary electrons, m e is the rest 
mass of the electron and c is the light speed. These primary 
electron-positrons will continue to emit curvature photons 
with a typical energy 



3 o c 
2^V 



6.4 x 10 8 P 1/2 eV, 



(2) 



where K is the Planck constant, s = (RRl) 1 / 2 is the cur- 
vature radius, R — 10 6 cm is the radius of the neutron 
star, and Rl = cP/2tt is the light cylinder distance. These 
curvature photons will be converted into secondary 
pairs in the local magnetic field when their energy sat- 
isfies (Ruderman & Sutherland 1975) 



E > E c 



2m e c 2 B q 
15 B~. 



3B~lGeV, 



(3) 



where B q = 4.4 x 10 13 G is a critical magnetic field, B St g 
is the local surface magnetic field of the neutron star in 
units of 10 9 G. 

It has been proposed that there is strong multipolc 
magnetic field near the stellar surface although a global 
dipole magnetic field gives a good description of the mag- 
netic field far from the star (Ruderman & Sutherland 1975; 
Ruderman 1991; Arons 1993). The typical radius of cur- 
vature I of the local magnetic field is of the order of the 
crust thickness of the star (i.e. I ~ 10 5 cm), which is much 
less than the dipole radius of curvature s of dipole field 
component near stellar surface. The relation between the 
local multipole magnetic field and dipole field can be given 
by (Zhang & Cheng 2003) 



B s ^B d (ff. 



(4) 



Bd is the dipole magnetic field of a pulsar. So the critical 
energy is rewritten as 



_ 2m e c 2 B qf R^_ 3 
crit " 15 B d [ l> ■ 



(5) 
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The secondary e ± pairs will lose their energy via syn- 
chrotron radiation with a photon energy at a distance of 
r 

E sy , = hl^f^m = 2 .6xlO«Ba, 9 P(**)»(£)»eV 
2 2m e c z m e c r I 

We can see that generally E syn > E CTlt , so a photon- 
electron cascade will start and develop until this condi- 
tion fails. At the end of a cascade, each incoming primary 
electron-positron can produce, on average 

N e± = -^ = 1.9 x 10 3 B d , 9 P 1 /3 ( ^ ) 3 ; (?) 

-frcrit I 

then the total pair production rate can be estimated as 
iV e± = fN Ne± = 2 x lO^/flyp- 8 / 21 ^) 30 ^- 1 , (8) 
where 

N = 2.7 x 10 27 P- 2 B lI ,9(y) 3 s- 1 (9) 

is the primary electron-positrons passing through the po- 
lar gap (Goldreich & Julian 1969). 

Since these pairs are created close to the stellar surface 
and the field lines arc converging, only a small fraction 
may keep moving toward the star and annihilate on the 
stellar surface. Ho (1986) showed that the loss cone for 
these pairs will approach ir/2 (cf. Eq. 29 of Ho 1986). 
In other words, most pairs will be reflected by magnetic 
mirroring effect and then move toward the light cylinder. 
These particles will flow out with the pulsar wind and be 
accelerated by the low-frequency electro-magnetic wave. 
In theory, the maximum Lorcntz factor of pulsar wind 
particles could be larger than 10 4 (Rees & Gunn 1974; 
Kennel & Coroniti 1984). 

3. Three possible scenarios for pulsar winds in GC 

As discussed in the introduction, we consider three pos- 
sibilities of positron sources from pulsar winds in the 
Galactic center region. In this section, we discuss them 
separately in detail. 

3.1. Normal pulsars 

The deep radio survey (LaRosa et al. 2000) and recent 
observations in the Galactic center field targeted in the 
deep Chandra X-ray surveys (Wang et al. 2002; Muno 
et al. 2003) showed that a few supernova remnants are 
present. Detailed analyses of X-ray nebulae GO. 13-0. 11, 
G359.89-0.08 and G359.54+0.18 (Wang, Lu & Lang 2002; 
Lu, Wang & Lang 2003) suggested young pulsars might 
lie inside them. So in the past years, young pulsars may 
have been produced actively in the Galactic center re- 
gion. However, because the average proper motion ve- 
locity of normal pulsars is high (~ 500 km s _1 ), few of 
them (< 40%) lie in the GC (Arzoumanian, Chernoff, & 
Cordes 2002), so the number of normal pulsars whose ac- 
tive timescale is ~ 1 Myr in the GC should be small, 



and the number of Crab or Vela like pulsars is around 10 
(Wang, Jiang & Cheng 2005). 

There is both theoretical and experimental evidence 
that young pulsars can emit intense electron/positron 
(^inds. In case of Crab and Vela, the injection rates 
are up to ~ 10 41 e + s" 1 and ~ 10 39 e + s _1 . Then in the 
timescale of the age of Crab or Vela, the total number 
of positrons is about 3 x 10 51 . Assume the birth rate of 
young pulsars is 1 per 10 3 years in the GC, and the life- 
time of positrons is 1 million years (cf. the cooling time 
scale estimate in section 4), there will be about 3 x 10 54 
positrons in the region. On average the annihilation rate 
is ~ 3 x 10 54 /lmillion yr ~ 10 41 s _1 . Therefore, our anal- 
ysis suggests that the young pulsar winds could not be 
a significant contributor to the positrons in the Galactic 
center region. 

3.2. Millisecond Magnetar in a GRB progenitor 

It has been suggested that gamma-ray bursts in the GC 
are potential positron sources (Casse et al. 2004; Parizot et 
al. 2005). Their models assumed that the annihilation line 
may relate to nucleosynthesis in the hypernova which is 
the GRB progenitor. But the GRB progenitor could also 
be a rapidly spinning strongly magnetized neutron star 
(millisecond magnetar, see Usov 1992; Thompson 1994), 
which is also a very strong pair emitter. The generic 
magnetized e ± outflow from a spinning compact object 
could be given by 

L m ~ 27rr 2 c(^-) ~ 1.5 x 10 52 r 2 6 B 2 5 (tA )- tt erg s" 1 , (10) 

ro is the size of the central engine, ~ 10 6 cm, B ~ 10 15 
G is the magnetic field of magnctars, t is a characteristic 
timescale after the burst, i.e. to ~ 1 sec, where we have 
assumed the magnetic field decay would follow a power 
law and a > 1. Assuming the average energy of positrons 
is ~ 1 GcV, the total number of positrons emitted by 
magnetars is 

roc 

N~ 2xl0 55 (t/t o y a dt~2xlQ 55 /(a-l). (11) 

J t 

If we take a = 1.1, the total number is up to 2 x 
10 56 . For other possible values of a, the total number 
will be reduced by a factor of several. On the other hand, 
the toroidal magnetic field of magnetars can be up 10 17 G, 
and part of this toroidal magnetic field energy can also be 
converted into e 1 *" (Kluzniak & Ruderman 1998; Dai & Lu 
1998). 

Therefore, in this model more than 10 56 positrons 
would be produced. Because the GRB birth rate from the 
present observations is ~ 10~ 5 yr -1 in a galaxy (Piran 
2004) after considering the beaming effect correction, we 
expect the GRB rate in the GC is about one per 10 6 years. 
Therefore, the annihilation rate of positrons produced by 
the magnetar is ~ 3 x 10 42 s _1 . But more comments on the 
GRB rates are required here. The observations of the long 
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GRBs and host galaxies are in agreement with a corre- 
lation of GRB locations with blue regions/spirals of host 
galaxies (Bloom et al. 2002), implying that long bursts do 
not concentrate in the bulge region; the theoretical stud- 
ies of merging neutron stars and black holes suggest that 
most short bursts may be located in the halo region (e.g 
Voss & Tauris 2003) , so here we possibly have overesti- 
mated the GRB rate in the Galactic center region consid- 
ering the uncertainties of burst sites of the different GRB 
populations. 



3.3. A millisecond pulsar population 

We have argued that millisecond pulsars could be the ma- 
jor contributor to the pulsar population in the Galactic 
center (see Wang, Jiang & Cheng 2005). There may exist 
a population of old neutron stars with low space velocities 
which have not escaped the Galactic center (see Belczynski 
& Taam 2004a; 2004b). Such neutron stars could be recy- 
cled to millisecond periods, and remain in the GC through 
their lifetime because of their relatively low proper mo- 
tion velocity (Lyne et al. 1998; Arzoumanian, Chernoff, 
& Cordes 2002). Wang , Jiang & Cheng (2005) suggest 
that the millisecond pulsar population can contribute sig- 
nificantly to the diffuse gamma-rays in the Galactic cen- 
ter (Mayer-Hasselwander et al. 1998). In addition, the 
X-ray synchrotron nebula around the millisecond pulsar 
B1957+20 (Stappers et al 2003) provides evidence that 
millisecond pulsars can remain active to emit relativistic 
winds. Therefore, we propose that electron-positrons from 
the winds of the millisecond pulsar population could be 
the possible positron sources in the GC. 

Since millisecond pulsars can lie in the GC and remain 
active for 1 Gyr, we believe they are continuous positron 
injecting sources. For the typical parameters of millisec- 
ond pulsars, P = 3 ms, fij ~ 3 x 10 8 G and accord- 
ing to equation 8, we estimate the positron injection rate 
N e ± ~ 5 x 10 37 s _1 for a millisecond pulsar. Wang et al. 
(2005) showed that about 6000 millisecond pulsars could 
contribute to the diffuse gamma-rays in the Galactic cen- 
ter region covered by EGRET. Since the present annihila- 
tion observations show the line emission region toward GC 
is about 6°, is much larger than the EGRET observation 
region (1.5°), so a larger number of MSPs are expected to 
lie in the bulge region of 6° . We do not know the distribu- 
tion of MSPs in GC, so we just scale the number of MSPs 
by 6000 x (6°/1.5°) 2 ~ 10 5 , where we assume the num- 
ber density of MSPs may be distributed as pmsp oc r~ x , 
where r c is the scaling size of GC. Based on the popu- 
lation analysis of Lyne et al. (1998), the number of mil- 
lisecond pulsars in the entire Galaxy may exceed 3 x 10 5 , 
so the number of millisecond pulsars in the bulge region 
estimated here could be reasonable. So the total positron 
injection rate from the millisecond pulsar population is 
- 5 x 10 42 c+ s" 1 . 



4. Electron-positron cooling in GC and 
annihilations 

The outflowing electron-positrons could be trapped by the 
magnetic field in the Galactic center. The Larmor radius 
{tl) of a relativistic electron with energy E e is given by 

r L « E e /eB ~ 10 13 (S e /10 10 cV)(B/10" 5 G)" 1 cm. (12) 

The angular size of annihilation line emission detected 
by SPI/INTEGRAL is several degrees (Knodlseder et al. 
2003), so we take A ~ 300 pc as the size of the Galactic 
center region we studied in this paper (we assume the 
distance of the GC to us ~ 8 kpc). The diffuse timescale 
of the electron is estimated as 

tdiff ~ (-) 2 - ~ 10 10 (£ e /10 10 eVrV, (13) 
r L c 

where we have taken the average magnetic field in the GC 
B ~ 10" 5 G (Uchida & Gusten 1995; LaRosa et al. 2005) 
and A ~ 10 21 cm. 

Electrons and positrons in the Galactic center will lose 
their energy via synchrotron radiation and Coulomb in- 
teractions with the medium ,and then annihilate into 511 
keV photons. If we assume the injection electron-positron 
energy spectrum as the following form: Q e ± ~ AE~ P , 
where p > 2, then positrons with energies higher than 
<~ 100 MeV will cool through synchrotron radiation in 
the magnetic field, while a large number of low energy 
positrons lower than ~ 100 MeV lose their energy by 
Coulomb interactions in the medium of the GC. The other 
energy loss processes, like inverse Compton cooling and 
brcmsstrahlung are negligible. When the positrons reach 
an energy around a few tens of eV, they could either an- 
nihilate directly with electrons or form positronium by 
charge exchange or radiative capture (Bussard, Ramaty 
& Drachman 1979). So the annihilation timescales depend 
on the cooling timescales. 

The synchrotron cooling timescale of relativistic 
positrons is estimated as 

r s - ~ 3 x lO^lO^eVr^V. (14) 

We find that the electron synchrotron timescale is lower 
than the diffusion timescale for the same electron energy, 
t s -C tdiff, so the positrons could significantly lose their 
energy to become non-relativistic, and stay in the Galactic 
center. 

The energy loss timescale of Coulomb interactions 
in the medium, t c depends on the positron energy, the 
medium density and the degree of ionization. For the 
positrons with energy around 1-100 MeV, the Coulomb 
interaction cooling timescale is estimated as r c <~ 10 5 n _1 
yr (within a factor of 2), and n is the medium number 
density in units of cm~ 3 . However, the medium in GC 
is quite complicated, there exist many giant molecular 
clouds (LaRosa et al. 2000), and hot gas discovered in X- 
rays (Muno et al. 2004). Then the energy loss timescales 
could vary from r c < 10 3 yr in the molecular clouds 
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(n > 10 2 cm~ 2 ), r c <~ 10 5 yr for the typical warm interstel- 
lar medium (n ~ 1cm -2 ), and t c > 10 7 yr for the hot gas 
(n ~ 10~ 2 cm -2 ). The analysis of the 511 keV annihilation 
spectrum suggests that the dominant fraction of positrons 
(~ 94%) form positronium before annihilation, which con- 
strains the positron annihilation in the hot gas to a very 
small fraction (< 8%, see Churazov et al. 2005). Most 
positron annihilations seem to lie in the warm medium, so 
the cooling timescale is also much lower than the diffusion 
timcscale. 



5. Summary and discussions 

The origin of Galactic positrons is still a mystery at 
present. In this paper, we have proposed a possible con- 
tribution to positron sources from the pulsar winds in the 
Galactic center region. We discussed the contributions of 
three possible pulsar scenarios in the GC: normal pul- 
sars; the rapidly spinning strongly magnetized neutron 
stars in GRB progenitors (millisecond magnetars); a mil- 
lisecond pulsar population. Electron-positron pairs are in- 
jected into the Galactic center region from three classes of 
pulsar winds, and could be trapped in this region by the 
magnetic field in the Galactic center. These relativistic e 
pairs will lose their energy though the synchrotron radi- 
ation and the Coulomb interactions with the medium to 
become non-relativistic. The cooling timescales are lower 
than the diffusion timescale of positrons in the magnetic 
field in the GC. So positrons could form positronium by 
charge exchange or radiative capture, or annihilate di- 
rectly with electrons into photons to produce the 511 keV 
line observed by the present detectors. Our results have 
shown that normal pulsars cannot be a significant con- 
tributor to positron sources in GC, but the millisecond 
magnctar and a millisecond pulsar population could be 
potential positron sources. However, as we have estimated 
that the diffusion time is much longer than the cooling 
time, it is not clear how the positrons produced by GRBs 
can spread over the entire galactic bulge or at least in the 
range of 6°. Parizot et al. (2005) have argued that the tur- 
bulent diffusion process is able to diffuse the positrons over 
the galactic bulge in a time scale of 10 7 yr. However, it is 
still not clear why the positrons produced by GRBs should 
concentrate in the galactic bulge, which is not supported 
by the observations of GRBs and their host galaxies (e.g. 
long bursts in Bloom et al. 2002; short bursts in Gehrels 
et al. 2005). In view of the uncertainty of the birth rate 
of GRBs in the GC and lacking of annihilation lines in 
disk, we have some reservations on GRBs as the signifi- 
cant positron sources in the GC. 

The predicted annihilation rate of positrons from the 
millisecond pulsar population is ~ 5 x 10 42 s _1 which 
is consistent with the present observational constraints 
by SPI/INTEGRAL. So we conclude that the e ± pairs 
from these classes of pulsar winds are potential positron 
sources, and can significantly contribute to the 511 keV 
annihilation lines in GC. 



The important problem for positrons is to determine 
the e + e~ annihilation line intensity and radiation mor- 
phology in the Galaxy. Recently, analyses of the data by 
SPI/INTEGRAL with deep Galactic center region expo- 
sure show that the spatial distribution of 511 keV line ap- 
pears centered on the Galactic center (bulge component), 
with no contribution from a disk component (Teegarden 
et al. 2005; Knodlseder et al. 2005; Churazov et al. 2005). 
We suggested that the positrons produced by pulsar winds 
could be trapped in the Galactic center by the magnetic 
field, the diffusion timescale is higher than the cooling 
timescales, so that the positrons can annihilate into 511 
keV photons before they can escape from the Galactic 
center region. In the other viewpoints, the bulge domi- 
nated morphology for 511 keV line may indicate that the 
positron source population could be an old stellar popu- 
lation. Knodlseder et al. (2005) have suggested that the 
low-mass X-ray binaries (LMXBs) could be the candidate 
sources because more than 60% of the Galactic LMXBs 
are observed towards the galactic bulge (see Grimm et 
al. 2002). The millisecond pulsar population is a very old 
stellar population, and part of them stay in LMXBs, but 
most may be isolated as radio pulsars, X-ray sources, even 
detectable in gamma-rays. The present radio observations 
find most millisecond pulsars lie in globular clusters (see 
the review by Camilo & Rasio 2005 and the references 
therein), few of them in the Galactic field. If assuming 
the millisecond pulsar evolutionary formation channels in 
globular clusters are similar to the Galactic center re- 
gion or bulge region, then we expect a millisecond pulsar 
population in the GC as a significant contribution to the 
positron sources. The millisecond pulsars can stay in the 
Galactic center region throughout their lifetime, which is 
well consistent with the 511 line emission morphology. 

We have shown that the millisecond pulsar popula- 
tion in the Galactic center could be the major sources 
of positrons. Normal pulsars could not be the main 
positron sources, but GRB progenitors, i.e. millisecond 
magnetars, may be significant contributors. In addition 
hypernovae/GRBs could also be the potential positron 
sources (Casse et al. 2004). Thus, how could we discrim- 
inate the model of a millisecond pulsar population from 
other models, specially models related to GRBs? Firstly, 
we can estimate the typical spatial diffusion scale, accord- 
ing to Eq. (13), Xdiff ~ (f'Lci) 1 / 2 , the average cooling time 
of positrons in GC is 10 6 years (Churazov et al. 2005), so 
the characteristic diffusion scale is about 10 18 cm. Because 
of the low angular resolution of SPI/INTEGRAL (about 
2 degrees, Vedrenne et al. 2003), we can assume that 
the positrons annihilate in the same local region as their 
sources, i.e. the millisecond pulsars. Therefore we predict 
that the spatial intensity distribution of the annihilation 
lines should follow the spatial distribution of MSPs if there 
exists a millisecond population in the GC. Additionally, 
we could assume the spatial distribution of MSPs should 
follow the mass distribution of the GC though we do not 
know how well they follow each other. But because the 
proper motion velocity of MSPs is relatively low, we could 
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reasonably assume two distributions are quite close to each 
other. Therefore, if the positron sources originate from the 
MSP population, the 511 keV annihilation line intensity 
would follow the mass distribution of the Galactic center 
region. The present image of 511 keV line emission (see 
Figs. 4 and 6 in Knodlseder et al. 2005) shows that the 
flux is strongest in the center, and decreases toward the 
outside region, possibly tracing the mass distribution of 
the GC, but this needs confirmation by high resolution 
observations in future. The millisecond pulsar population 
could naturally explain the diffuse emission morphology 
and does not have the problem of the turbulent diffusion, 
which is required to diffuse all these positrons to a few hun- 
dred pc. On the other hand, for models associated with 
GRBs, the positron source is initially a point-like source, 
so the diffuse emission of 511 keV line needs the strong tur- 
bulent diffusion to diffuse the positrons to a larger scale 
(Parizot et al. 2005). In the present viewpoints, GRB pro- 
genitors (either magnetars or hypernovae) are related to 
massive stars and dense clouds (e.g. molecular clouds in 
GC which has a dimension of a few tens pc, LaRosa et al. 
2000), the mean density n > 10 2 cm~ 3 , the cooling time 
of positrons is less than 10 3 years. So the spatial diffu- 
sion scale is only 10 17 cm, which makes the annihilation 
line emission looks like a point source correlated with the 
molecular clouds. In order to make the line emission be a 
diffuse source, GRB positrons must diffuse to a hundred 
pc with a long diffusion timescale > 10 years, so it re- 
quires that the mean density in the GC is ~ 10~ 2 cm~ 3 . 
But this low density seems unreasonable in the GC, and 
also is not consistent with the constraint on the medium 
environment from the spectrum of 511 keV line (Churazov 
et al. 2005). Furthermore, if positrons cannot diffuse to 
larger distances, the intensity distribution of annihilation 
lines should not correlate with the general mass distribu- 
tion of the GC, instead it should be more correlated to the 
distribution of massive stars and dense molecular cloud. 
Moreover, GRBs more likely occur in the disk, where more 
dense molecular clouds and massive stars are located, and 
in fact the detected magnetars and soft gamma-ray re- 
peaters are also distributed in the disk. Giant flares from 
soft gamma-ray repeaters (SGR) are also positron source 
candidates which are distributed in the Galactic plane, but 
the 511 keV line observations show a very weak disk com- 
ponent (Teegarden et al. 2005; Knodlseder et al. 2005). 
The 27 Dec. giant flare from SGR 1806-20 released an en- 
ergy up to 10 47 erg (Hurley et al. 2005), and is thought 
to be a possible origin of short GRBs, and could be ob- 
served as the 511 keV emission source in future. In sum- 
mary, the positron source models associated with GRBs 
are disputable. Our scenario of a millisecond pulsar pop- 
ulation as possible positron sources in the GC has some 
advantages to explain the diffuse morphology of 511 keV 
line emissions, and predicts the line intensity distribution 
should follow the mass distribution of the GC, which may 
be tested by future high resolution observations. 
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